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A NOTE ON THE ABILITY TO PREDICT TRANSONIC : 
CHANGES DUE TO MODEL MODIFICATIONS 
By P. Kenneth Pierpont 


SUMMARY 


An experimental and 'analytical study was made at transonic speeds up 
to a Mach number of 1.43 "to ascertain the ability of present calculation 
procedures to predict transonic drag-rise changes which result from phys- 
ical model changes. The experimental data were obtained in the Langley 
8-foot transonic tunnel and the Langley 8-foot transonic pressure tunnel 
at Reynolds numbers based on the mean aerodynamic chord near 2 X 10^. 

Both wing-body and wing-body-tail configurations were investigated. 


A 

)RAG-RISeQ 



The results showed that, with the exception of a model having a 
relatively "bluff forebody and afterbody, the use of current techniques 
utilizing a Fourier analysis to fit the area slope curve of the equiva- 
lent body provided estimates of transonic drag within 15 percent of the 
measured values for wing-body-tail configurations. When a well-shaped 
configuration was chosen as a basis, present calculative procedures per- 
mitted estimates of the change in drag rise which resulted from physical 
change in the model which were considered accurate enough to be useful 
for preliminary design. 


INTRODUCTION 


The application of the transonic and supersonic area rules (refs. 1 
and 2) aroused considerable interest in methods for numerically estimating 
transonic and supersonic wave drags. Computing techniques such as those 
outlined in references 3 and 4 have been effectively used for estimating 
the zero-lift drag-rise coefficients at transonic speeds of model config- 
urations (refs. 5 and 6, for example) and for complete configurations in 
reference 7 . 

The ability to predict transonic drag-rise coefficients above M » 1.02 
closer than the 20-percent value stated In reference 3 depends in part at 
least upon the experience of the computer and the characteristics of the 



2 


NACA RM L57H07 


individual configuration. For those cases in which performance estimates 
demand more accurate information, resort to wind-tunnel or free-flight 
model tests ."becomes a necessity. However, since a given "basic configura- 
tion may change considerably as the design progresses, some method for 
estimating the effects of external modifications is needed. The purpose 
of the present study was to determine if the changes in drag rise asso- 
ciated with physical modifications can be adequately predicted by the 
numerical procedures presently available. 

A swept wing and tail model, typical of current fighter designs, was 
tested in the Langley 8-foot transonic tunnel and in the Langley 8-foot 
transonic pressure tunnel. Various arrangements of fuselage modifications 
together with wings and tails of different thicknesses were employed. The 
Mach number range extended from subcritical speeds to 1.2 and in some 
cases to 1.4, and the Reynolds number based on the wing mean aerodynamic 
chord was about 2 x 10^. 


SYMBOLS 


A 

An 

CD 

ACp 

ac d > 

l 

M 

N 

q. 

s 

X 


cross-sectional area, normal to free-stream direction 

Fourier coefficient of nth term 

drag coefficient, Drag/qS 

drag-rise coefficient, - Cp^^o.8) 

incremental drag-rise coefficient, corresponding to drag rise 
of an equivalent body having an area distribution for a 
particular Mach number at a particular value of 0 

body length 

Mach number 

number of terms used in Fourier sine series 
dynamic pressure 
wing area 

distance along X-axis, measured from body nose 


■ * 
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angle between Z-axis and line in YZ-plane normal to intersection 
of cutting plane with YZ-plane (see fig. l) 

_ / Oy \ 

longitudinal position angle, cos (1 - — J 

MODEL AND TESTS 


The model consisted of a nonaxi symmetric fuselage having the wing 
mounted on top with 42° of quarter-chord sweepback and an aspect ratio 
of 3.4. Model details are shown in figure 2 and table I. Tail surfaces, 
when used, had 45° quarter-chord sweep and aspect ratios of 2.0 and 1.5 
for the horizontal and vertical tails, respectively. 

Two wings, made up of straight-line elements, were employed; the one 
having a taper in thickness from 6 to 5 percent was called the basic wing 
and the other of a uniform 4-percent thickness was called the thin wing. 
Similarly two sets of tail surfaces were used, the basic tails having a 
taper in thickness from 6 to 4 percent and the thin tails having a uni- 
form 4-percent thickness. Symmetrical airfoil sections of the NACA 
65 A series were used for all surfaces. 

Several fuselage configurations, designated by letters, are shown 
in figure 5 . A configuration described by letter only signifies a basic 
wing and body; whereas, a letter with a prime signifies a body with basic 
wing and tails, and the subscript 1 with either the letter or the letter 
and prime indicates that the 4 -percent -thick surfaces have been installed 
in place of the tapered-in- thickness surfaces. For example: C designates 

body C with basic wing, C 1 designates body C with basic wing and basic 
tails, and C'i designates body C with thin wing and thin tails. 

The experimental investigation was conducted in both the Langley 
8 -foot transonic and the Langley 8 -foot transonic pressure tunnels 
described in references 8 and 9* A sketch of the model in the Langley 
8 -foot transonic pressure tunnel is shown in figure 4. Mach numbers 
extended from 0.8 to about 1.2 with a few points being obtained at 1.43. 
Reynolds number based on the wing mean aerodynamic chord was about 

2.0 X 10^. 

Zero-lift drag was obtained from faired polars for a small range of 
angles near zero lift. Estimated accuracy for the drag-coefficient data 
is considered to be about ±0.001. Although there is some difference in 
the turbulence characteristics of the two tunnels, unpublished data indi- 
cate that the difference in zero-lift drag coefficient measured in both 
tunnels is approximately constant throughout the range of test Mhch num- 
bers. Comparisons of drag-rise coefficients obtained in either facility 
should therefore be valid. 
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COMPUTATIONS 


The computations necessary to obtain the drag-rise coefficients were 
carried out in the manner outlined in reference 3* The drag-rise coef- 
ficient for a complete configuration was defined as 


ac d = 



( 1 ) 


in which AC^ ' is termed the incremental drag-rise coefficient and is the 
calculated drag-rise coefficient for the equivalent body corresponding to 
a particular single area distribution at the angle <jf. (See fig. 1.) 
Accordingly, from reference 3> 


'D 


-r s l< 


n=l 


( 2 ) 


Values of A n were obtained by conventional Fourier analysis techniques 
on a digital computer and are defined as 


An 



sin n9 d9 


(3) 


The wing and tail empennage area distributions were obtained by 
graphical means. This method was found to be both time consuming and 
subject to appreciable human inaccuracies. Exploration of numerical 
methods for obtaining these areas led to solutions which were basically 
similar to those in references 10 and 11. Check computations by the 
numerical methods indicated that satisfactory results could readily be 
obtained and therefore served as a useful check on the graphical solutions. 

Area distributions normal to the longitudinal axis of the fuselage 
were obtained by integrating photographically reproduced cross sections 
cut from plaster of Paris mold patterns. In the region of the wing-body 
juncture as well as in the region of the tail-body junctures a small 
portion of the surface was arbitrarily included in the defined fuselage 
area to simplify the work. 
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Total equivalent-body area distributions were obtained by combining 
the various components as outlined in reference 12. In this reference, 
it was suggested that the cross-sectional area for the wing or horizontal 
tail be combined with the fuselage areas at the longitudinal stations 
where the oblique cutting planes for the wing or tail cross the plane of 
lateral symmetry. Cross-sectional areas for the vertical tail were com- 
bined with the fuselage areas at the stations where the cutting planes 
crossed the upper surface of the fuselage near the base of the tail. 

Slopes of the area distributions were derived numerically by using 
a five-point analysis for 120 points equally spaced along the body length. 
A check on the validity of the slope curves was obtained by plotting and 
integrating mechanically to insure the proper closed body. Thirty-three 
harmonics were computed for all configurations after it had been noted 
that the area distribution was not adequately represented by the initially 
selected 2k harmonics. 


‘.RESULTS AND DISCUSSION 


The availability of data for several body configurations tested in 
conjunction with different wing and tail thicknesses has enabled an evalu- 
ation of available methods for computing the drag-rise -coefficient changes 
which resulted from diverse model modifications. For convenience of 
presentation, the results have been arranged in two categories. The first 
group, figures 5 and 6 , consists of wing-body configurations and includes 
some of the effects of wing thickness. The second group, figures 7 and 
consists of wing-body-tail configurations and includes additional effects 
of wing and tail thickness. Representative variation of the incremental 
drag-rise coefficient, which is the drag rise for a body having the area 
distribution equivalent to that for a particular cutting-plane angle 0 , 
with the cutting-plane angle is shown on figure 9 . The tapered-in- thickness 
wing and tail will hereinafter be referred to as the basic wing and tails 
and the uniform thickness 4-percent wing and tails will be termed the thin 
wing and tails. 


Wing-Body Configurations 

Body shape .- Longitudinal equivalent -body area distributions for 
three wing-body configurations with the basic wing and two wing-body 
configurations with the thin wing are shown in figure 5 for three prin- 
cipal cutting planes. Model A is shown to have a blunt forebody shape 
and a large negative slope near the base. The representative area slope- 
distribution curve shown in figure 6 illustrates this more readily. 
Substantial reductions in local slopes at both nose and tail were achieved 
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by lengthening and slenderizing the forebody and by adding a small exten- 
sion to the model base; these changes are illustrated by body B (see also 
fig. 3). 

Further improvement in the overall area distribution was achieved by 
filling in the hollow ahead and behind the maximum area of body B to 
obtain model C. These area additions were obtained from an approximate 
average area distribution for M = 1.2 by using a tangent straight line 
across the hollow as described in reference 11. The area to be added was 
then divided so that approximately one-half was added above and one-hali* 
below the wing chord plane. The representative slope-distribution curves 
of figure 6 show that this modification has reduced the magnitude of the 
peak slopes. Also shown in the figure are a number of check points used 
to establish the degree of approximation for body B obtained from the 
Fourier analysis. These check points are representative of the agreement 
achieved throughout the analysis. 

Comparisons of calculated and experimental drag-rise coefficients 
are shown in figure 10(a) . It appears plausible that the lack of agreement 
for model A could be attributed to the inability of the theory to properly 
account for the relatively bluff forward and rearward portions of the 
body. The agreement between calculated and experimental results for 
models B and C are excellent. The estimated change in drag-rise coeffi- 
cient between models B and C agreed within about 30 percent with the 
measured values at both M = 1.0 and 1.2. 

Wing thickness .- The two models, Bi and Fp, in figure 5 were 
equipped with the thin wing of 4-percent thickness. Body F was nearly 
identical to body C. Principal effect of reducing the wing thickness 
on the equivalent -body geometry was to increase the equivalent-body fine- 
ness ratio which can be seen from figure 5« Because body F had been 
designed with the basic wing, some overfilling of the hollows ahead and 
behind the maximum area is evident when used with the thin wing, especially 
at <f = 0° and 90°. For these configurations only the calculated drag- 
rise coefficients were obtained. Values of C D = 0.0146 and 0.0133 
correspond to configuration Bp and Fp, respectively, for M = 1.2. 
Compared to configuration B, the change in drag-rise coefficient resulting 
from the wing modification (model Bp) would be 0.0071 and compared to C 
the change resulting from the wing modification (model Fp) would be 
only 0.0024. 


Wing-Body-Tail Configurations 

Body shape .- Longitudinal equivalent-body area distributions for 
five wing-body-tail configurations using the basic wings and tails are 
shown in figure 7 . Two configurations for which the thin wings and tails 
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were employed are compared in figure 8. In the latter figure, body F 
was nearly identical to body C. 

Comparison of the calculated values of incremental drag coefficient 
as a function of circumferential position angle or cutting plane angle 
are shown in figure 9 for configurations A f , B 1 , and C*. Variation with 
0 is appreciably different for the three models. Model C especially 
bears little resemblance to either of the others. Different variations 
with (f> such as are shown were typical of many of the configurations. 
Caution should therefore be exercised when attempting to predict overall 
drag-rise changes on the basis of area distributions for a single cutting- 
plane angle. 

Experimental drag-rise coefficients have been compared with the 
corresponding calculated values in figure 10(b) for both M = 1.2 and 1.37. 
As was the case for the wing-body configuration alone, model A 1 is shown 
to be overestimated by a significant amount. The calculations were in 
substantial agreement with the test data for models B 1 , C ! , D r , and E*. 
Model D* was obtained by deliberately overfilling the area distribution 
to obtain a shape, based on the average area curve, approximating an 
ideal-slender-body-theory shape as outlined in reference 11 near the 
maximum area. Model E was obtained by first indenting the fuselage under- 
neath the wing only in the region of the maximum area by about 10 percent 
of the equivalent -body maximum area, an amount considered a maximum with- 
out interfering with an engine, and then filling in the remaining depres- 
sion in the area distribution curve in a manner similar to that for con- 
figuration CM. Inconsequential improvement for model E 1 as compared to 
model C 1 was measured and the calculations actually showed a small drag 
increase. The calculated drag rise for model E 1 was the only actual 
inconsistency in the calculations and is not surprising since the theory 
cannot be expected to account for very small drag reducing changes in 
configurations. A similar case of asymmetric indentation was shown to 
have an adverse effect on the drag rise for a delta wing in reference 13. 
The decrease in friction drag with Mach number, however, would tend to 
improve the agreement between computed and experimental results particu- 
larly at the higher Mach number (M = 1.37) • 

The effects of body-shape change on the calculated drag rise with 
the thin wings and tails (fig. 8) are shown in figure 10(c) compared to 
that for model B 1 . For configurations B 4 ^ and F*-j_, the body shapes were 
the same as those previously used with the thick wings and tails. These 
shapes, therefore, do not necessarily represent near-optimum shapes but 
still are useful for comparison. Figure 10(c) indicates that the calcu- 
lated and measured values are in fair agreement for model B 1 j_ and in good 
agreement for model F'q. Furthermore, the improvement which would be 
predicted did not materialize. 
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For the design of an actual airplane, total drag instead of drag 
rise is the important parameter. For the present wing-body- tail config- 
urations, an average subcritical drag coefficient was found to he 0 . 015 - 
By using this value together with the calculated drag-rise coefficient, 
the error in estimated transonic drag was less than 15 percent except 
when body A was utilized. 

Wing and tail thickness .- Generally, the effect of reducing the wing 
and tail thicknesses was to reduce the component drag-rise contributions 
and to increase the overall fineness ratio of the equivalent body. Com- 
parison of figure 10(c) shows that, for body B, the reduction in wing 
and tail thickness resulted in drag-coefficient decreases which were 
adequately predicted to be about 0.010. Comparing model C' and F'q for 
which bodies C and F were nearly identical, shows that the calculated drag- 
rise coefficients agreed well with the measured values. The change in 
drag coefficient due to thickness change was only about one-half as much 
for body F as with body B, or about O.OO 5 . This comparison therefore 
indicates that the improved body shape had eliminated much of the inter- 
ference drag and that the remaining change was largely a result of the 
wing and tail wave drag. 


Drag-Rise Change 

The results have generally shown that the calculated and experimental 
results are in better agreement for those models having improved area 
distributions. The results have been replotted in figure 11 in the form 
of drag-rise changes. Drag-rise change is defined as the increase or 
decrease in drag rise resulting from a configuration change whether 
measured or calculated. For example, the drag-rise change for model C' 
referred to the base of B' is j B * ~ Cd,C'| . Both the tail-off and 
tail-on configurations have been compared on the basis of model B, whereas 
tail-on configurations have been compared on the basis of model B'. The 
amount of overestimation or underestimation resulting from the calcula- 
tions is shown by appropriate shading, and, except for the model with 
body A, the method is considered accurate enough to be useful for pre- 
liminary engineering design studies. 


CONCLUSIONS 


An experimental and analytical study has been made in the Langley 
8-foot transonic tunnel and in the Langley 8-foot transonic pressure 
tunnel to ascertain the ability of present calculation techniques to 
predict transonic drag-rise changes. The following conclusions are 
believed applicable: 
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1. With the exception of a model having a relatively "bluff forebody 
and afterbody, the use of currently available calculative procedures pro- 
vided estimates of transonic drag for wing-body-tail configurations 
within 15 percent of the measured values. 

2. When a well-shaped configuration was chosen at the outset, present 
calculative procedures permitted predictions of the change in drag rise, 
resulting from physical change in the model, which are considered accurate 
enough to be useful for preliminary design. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., July 18, 1957* 
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TABLE I.- MODEL DETAILS 

Wings: 

Sweepback, quarter-chord, deg 42 

Taper ratio 0.25 

Aspect ratio 3.4 

Mean aerodynamic chord, in , t 5.95 

Incidence at root, deg -1 

Twist, deg 0 

Airfoil sections: 

Basic root NACA 65 AOO 6 

Basic tip NACA 65 AOO 5 

Thin root NACA 65A004 

Thin tip NACA 65 AOO 4 

Horizontal tails: 

Sweepback, quarter-chord, deg 45 

Taper ratio 0.15 

Aspect ratio, based on total area 1.97 

Airfoil sections: 

Basic root NACA 65 AOO 6 

Basic tip NACA 65 AOO 4 

Thin root NACA 65 AOO 4 

Thin tip NACA 65A004 

Vertical tails: 

Sweepback, quarter -chord, deg 45 

Taper ratio 0.26 

Aspect ratio 1.47 

Airfoil sections: 

Basic root NACA 65 AOO 6 

Basic tip NACA 65A004 

Thin root NACA 65A004 

Thin tip NACA 65A004 

Fuselages : 

Length (without fairing) of configuration A, in 25.30 

Length of configuration B, in 27.88 

Fineness ratio of configuration A 8.63 

Fineness ratio of configuration B 9*^5 
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Figure 7*~ Comparison of the total-area distributions for the complete model. M = 1.2. 

(All dimensions are full scale.) 
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Figure 8.- Effect of empennage thickness on the area distributions. M = 1.2. (All dimensions 

are full scale.) 
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Figure 9.- Variation of incremental drag-rise coefficient with circum- 
ferential position. M = 1.2. 
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(a) Tail off; basic wing. 


Figure 10.- Comparison of the calculated and experimental drag-rise coefficients with Mach 

number . 











NACA - Langley Field, Va. 


VJ4 

O 


Configuration 
A C 

. 020 “! 


.015 



Base: B 

(a) Tails off. 


Configuration 



Base: B 1 



Overestimate 

Underestimate 


(t>) Tails on. 



Figure 11.- Comparison of estimated drag-rise increments. M = 1.20. 





